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Steady -State  Multi- Ion  Disturbed  D- Region  Model 


I.  I INTRODUCTION 

P eviousl v 1 we  described  a model  for  the  nighttime  disturbed  D-region  which 

2 

treated  negative  ions  in  a highly  simplified  manner.  Good  agreement  was  attained 
between  the  model  and  the  nighttime  observations  of  3 November  1969  at  Ft.  Churchill, 
Canada,  during  a solar  proton  event  (SPE).  More  recently,  by  incorporating  the 

3 

latest  negative  ion  chemistry,  exclusive  of  hydration  processes,  we  have  found'  that 
good  agreement  can  be  obtained  between  our  latest  model  and  the  three  daytime 
observations  of  2,  3,  and  4 November  1969  for  this  same  event  and  location.  In  this 
report  we  document  the  new  program  which  includes  not  only  a more  satisfactory 
treatment  of  negative  ions,  but  also  includes  some  new  positive  ion  chemistry  involv- 
ing NOf  ions.  We  also  provide  here  the  numerical  results  obtained  in  applying  this 
program  to  the  solar  proton  events  of  2 to  5 November  1969  and  3 to  1 1 August  1972. 


(Received  for  publication  14  June  1978) 

1.  Swider,  W.  , and  Foley,  C.  I,  (197  5)  Computer  Program  for  the  Disturbed 

Steady -State  Nighttime  D-region,  A FCRL  -Tft  -V  5 -0  l5o. 

2.  Swider,  W.  , and  Narcisi,  R.  S.  (1975)  A study  of  the  nighttime  D-region  during 

a PC  A event,  J.  Oeophys.  Res.  80:655-664. 

3.  Swider,  W.  , Keneshea,  T.J.,  and  Foley,  C.I.  ( 197  8)  An  SPE -disturbed 

D-region  model.  Planet.  Space  Sci.  26: 
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2.  OUTLINE  OK  MODEL 


The  computer  program  is  reproduced  in  Appendix  A.  All  photorates  in 
this  code  are  multiplied  by  DD  which  is  set  to  1 for  daytime  and  0 for  night- 
time. 

The  iterative  scheme  developed  to  solve  for  all  species  concentrations  is  simi- 
lar to  that  derived  before.  The  present  scheme  is  more  complex;  two  initial  con- 
centrations are  now  needed,  one  for  electrons  and  one  for  the  positive  ion  sum. 

Using  the  initial  concentrations  for  electrons  [e],  and  the  positive  ion  sum,  SP,  the 
individual  negative  ions  are- determined  and  then  summed,  the  sum  denoted  by  NSUM. 
Then  the  positive  ions  are  individually  calculated  and  summed,  PSUM.  A new 
[e]  = Jprevious  [e]  + PSUM/(1+A)|  / 2,  where  A = NSUM/previous  [e],  and  a new 
SP  = Jprevious  [e]  + NSUM  + PSUM  j / 2 are  then  determined.  This  sequence  is 
repeated  until  PSUM  agrees  with  NSUM+[e]  to  some  desired  precision,  usually 
± 1 percent.  In  the  examples  to  be  shown  later,  not  more  than  9 iterations  were 
required  at  any  particular  altitude  (40-90  km)  in  order  to  attain  this  criterion. 

Initial  daytime  electron  concentrations  were  obtained  from 

[e]  = (q/Sl)1/2  (1)  L 189 

where  SI  = > p represents  values  for  the  effective  electron  loss  coefficient  as  deter- 

4 

mined  for  the  2 to  5 November  1969  SPE.  The  number  following  the  equation 
number  in  the  text,  as  L 189  above,  means  that  Eq.  (1)  is  or  begins  on  line  189  of 
the  program.  Initial  nighttime  electron  concentrations  were  estimated,  as  before,  * 
from  the  equation 

(e J = J(L(A)/2aD>  +q/aDJ  -L(A)/2»D  (2)  L 191 

where  « ^ is  the  mean  dissociative  recombination  coefficient  and 

L(A)  = kgl|02]2  + kg2[02]  [N2]  (3)  L 187 

is  the  loss  rate  (sec  S of  electrons  by  attachment  to  02  with  O 2 and  Ng  as  third 
bodies.  The  reaction  rates  and  rate  coefficients  are  listed  later,  the  subscripts 
for  k referring  to  the  processes  are  as  numbered  in  that  listing 

4.  Swider,  W.  , and  Dean,  W.  A.  (1975)  Effective  electron  loss  coefficient  of  the 

disturbed  daytime  D-region,  J.  Geophys.  Res.  80: 1815-1819. 

5.  Swider,  W.  , Narcisi,  R.S.,  Keneshea,  T.  J. , and  Ulwick,  J.  C.  (197  1)  Electron 

loss  during  a nighttime  PCA  event,  J.  Geophys.  Res.  76:4691-4694. 
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Initial  total  positive  ion  concentrations  were  determined  from  the  expression 


SP 


/ 5iOJ_ 
\®D 


l°3]  \ / 

^r)/(5 


1/2 


5(0]  + [Og]) 


(4)  LI 93 


where  a.  is  the  mean  ion-ion  recombination  coefficient.  The  multiplicative  factor 
for  [O],  somewhat  arbitrary,  was  selected  after  some  trial  runs  as  a factor  for 
improving  SP  vis-a-vis  the  final  calculated  positive  ion  sum.  liquation  (4)  simplifies 

to  the  familiar  expressions  (q /or.  ) 1 ^ “ at  high  l)-region  altitudes  where  (0]>>[0,] 

1/2  J J 
and  (q/ar.)  at  low  D-region  altitudes  where  ([0„]>>[0],  Values  chosen  for  or. 

1 -73  ”8*^3  ^ 

and  Of.  were  4X  10  cm  /sec  (L  73)  and  6 X 10  cm  /sec  (L  74)  respectively.  The 

former  value  was  selected  as  that  typical  for  a mixture  of  NO+  and  ions  at 

D-region  temperatures,  whereas  the  latter  value  is  based  upon  laboratory  measure- 

g 

ments  by  Smith  et  al.  The  fact  that  all  ion- ion  recombination  coefficients  were 

“83  * 

set  at  6 X 10  cm  /sec  in  this  program  is  in  accord  with  the  laboratory  work.  How- 
ever, this  constraint  is  not  a necessary  one  in  the  operation  of  this  program;  one 
n?.  being  required  only  for  initial  SP. 

The  positive  ion  equations  are  similar  to  those  listed  before  except  that  changes 

have  been  made  because  of  our  inclusion  of  photodissociation  processes  for  0^+  and 

02+  (HgO)  ions  and  by  our  addition  of  the  ions  NO+  (N2>,  NO+(H20)  (N9),  NO+  (H20)(C02), 

NO+  (2H20)(N2>  and  NO+  (211^0)  (CO^).  The  list  of  positive  ion  reactions,  numbered 

4 through  52  in  Table  1 for  the  21  positive  ions  considered,  is  fairly  complete.  Rate 

2 3 

coefficients  are  generally  from  the  same  sources  already  identified.  ’ Rates  for 

7 

photodissociation  processes  4 and  5 in  Table  1 are  from  Vanderhoff  and  Niles, 
while  the  rate  coefficient  k^  was'taken  from  the  work  of  Johnsen  et  al.  ^ Rate 
coefficients  for  reactions  45  and  48  were  set  equal  to  that  for  k4g.  The  rate  coeffi- 
cient k51  is  based^  upon  parameters  measured  by  Johnsen  et  al.  8 High  rate  coeffi- 
cients for  processes  44,  46,  47,  49  and  50  were  chosen  in  accordance  with  the 
fact  that  they  are  all  "switching"  reactions.  Computation  of  the  individual  positive 

ions  begins  with  09+  (PION  4,  J),  line  204  of  the  program  and  is  completed  by 
^ + 

line  257  with  PION  (24,  J)=HgO^  . Thus,  this  rather  extensive  list  of  positive 
ions  (21)  requires  only  53  lines  in  the  computer  program. 


6.  Smith,  D.  , Adams,  N.G.,  and  Church,  M.J.  (1976)  Mutual  neutralization  rates 

of  ionospherically  important  ions.  Planet.  Space  Sci.  24:697-  703. 

7.  Vanderhoff,  J.  A.  , and  Niles,  F.  E.  (1976)  Recent  development  in  measurements 

of  photodissociation  for  atmospheric  positive  ions,  EOS,  57:303. 

vw 

8.  Johnsen,  R.  , Huang,  C.-M.,  and  Biondi,  M.A.  ( 1975)  The  formation  and 

breakup  of  NO+,N9  clusters  in  N9  at  low  temperatures,  J.  Chem  Phys. 

63:3374-3378.  Z 

w 

9.  Heimerl,  J.  M.  (1976)  Values  for  Selected  Collisional  Dissociation  Rate 

Coefficients  of  Positive  Clustered  Ions,  BRL-1MR-499. 
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Table  1.  Positive  Ion  Reactions.  Units  for  photo-processes,  two-body 
and  three-body  processes  are  sec'1,  cm3/sec  and  cm6/sec,  respectively. 
The  various  process  coefficients,  <c  = A (T/300)®  e D/T,  are  presented  in 
the  format  A B D 


A B 


I* 

04* 

♦ 

HJ 

s 

02* 

♦ 

02 

3.00E-01 

5 

02*.H20 

♦ 

Nyf 

= 

02* 

♦ 

H20 

6.00E-01 

S 

02* 

♦ 

HO 

= 

NO* 

♦ 

02 

4.50E-10 

7 

02* 

♦ 

02 

♦ 

02 

S 

04* 

♦ 

02 

3.90E-30  -3.2 

9 

0 4* 

♦ 

0 

= 

02* 

4- 

03 

3.00E-16 

9 

04* 

4 

H20 

= 

02*H20 

4- 

02 

1. 50E-09 

10 

0 2*H23 

4 

H20 

H30*5H 

♦ 

02 

1. DOE- 09 

11 

02*H20 

4 

H20 

= 

H30* 

♦ 

OH 

4 

02 

2.00E-10 

12 

H 30  + 0H 

♦ 

H 20 

= 

H502* 

4 

OH 

1.40E-09 

13 

H30* 

♦ 

H2 

f 

N 2 

= 

H30*.N2 

♦ 

N2 

1.40E-30  -2.0 

1«* 

H33+.N2 

4 

CO  2 

= 

H30* . C02 

♦ 

N2 

1.00E-09 

15 

H33* .C02 

4 

H20 

= 

H502* 

♦ 

C02 

1.00E-10 

IS 

H502* 

♦ 

C02 

♦ 

N 

= 

H502* • C02* 

N 

3.  OOE-30 

12 

H592+.C02* 

H20 

= 

H703* 

♦ 

C02 

1.0CE-10 

13 

H 30* 

4 

H20 

4 

i 

= 

H502* 

♦ 

M 

3.50E-27  -3.0 

19 

H502* 

4 

H 20 

4 

N 

s 

H703* 

♦ 

H 

2.20E-27  -3.0 

20 

H703* 

4 

H 2 0 

f 

M 

= 

H904* 

♦ 

H 

2.30E-27  -3.0 

21 

NO* 

♦ 

CO  2 

♦ 

N2 

= 

N0*C0  2 

♦ 

N2 

3.  00E-29  -2.0 

22 

N 0*C02 

♦ 

N20 

= 

N0*H20 

♦ 

C02 

5.  CQE-10 

23 

N0*H20 

♦ 

H20 

4 

N2 

= 

NO* ( H 20) 2* 

N2 

1.10E-27  -4.7 

24 

NO* (H20) 2 ♦ 

H20 

4 

N2 

= 

N0MH20)  3* 

N2 

1.60E-27  -4.  7 

25 

NO* (H  20) 3* 

H 20 

H703* 

4 

HN02 

7.00E-11 

26 

NO* 

♦ 

E 

= 

N 

4 

0 

4.00E-07  -1.0 

27 

0 2* 

V 

E 

= 

0 

4 

0 

2.00E-07  -1.0 

28 

0 4* 

♦ 

E 

= 

02 

4 

02 

2.00E-06 

29 

0 2* .H20 

♦ 

E 

= 

02 

4 

H20 

2.00E-06 

30 

H30*.0H 

♦ 

E 

= 

2H20 

2.00E-06 

31 

N 30* 

♦ 

E 

H2 

4 

OH 

1.00E-06 

32 

H 502* 

♦ 

E 

= 

H20 

4 

H2 

4 

OH 

2. 00E-06 

33 

H703* 

♦ 

E 

= 

2H20 

4 

H2 

4 

OH 

4.00E-06 

34 

H904* 

4 

E 

= 

3H20 

4 

H2 

4 

OH 

5. 00 E- 06 

35 

H30*.N2 

♦ 

E 

3T 

H 

4 

H20 

4 

N2 

2.00E-06 

36 

H30*. C02 

♦ 

c 

= 

H 

4 

HZO 

4 

C 02 

2. 00E-Q6 

37 

H502* .C02* 

E 

= 

H 

4 

2H20* 

C 02 

3.2  0E- 06 

38 

N 0*  • C 0 2 

♦ 

E 

= 

NO 

4 

C02 

1.50E-06 

39 

N 0*H20 

4 

E 

= 

NO 

4 

HZO 

1.50E-06 

40 

NO*  (H20)  2 ♦ 

E 

S 

NO 

4 

ZH23 

3.00E-06 

41 

NO* (H20I 3* 

E 

= 

NO 

4 

3H20 

4.50E-06 

42 

04* 

♦ 

H 

= 

02* 

4 

02 

♦ N 

3.30E-06  -4. 0-5.03E*  03 

43 

NO* 

♦ 

H 2 

f 

N2 

NO* • N 2 

4 

N2 

2.00E-31  -4.4 

44 

N0*.N2 

♦ 

CO  2 

= 

NO* . C 02 

4 

N2 

2.00E-09  0.0 

45 

N0*.H20 

♦ 

N2 

♦ 

N 2 

= 

N0*H2ON2 

4 

N2 

2.  0 0E-31  -4.4 

46 

N 0*H20N2 

♦ 

C 02 

= 

N0*H20C02* 

N2 

2.00E-09  0.0 

47 

N0*H20C02* 

H20 

NO  * . 2 H20 

4 

C02 

2.00E-09  0.0 

48 

NO*. 2H20 

♦ 

N 2 

4 

N 2 

= 

NC*2H20N2* 

N2 

2.00E-31  -4.4 

49 

N0*2H20N2* 

C02 

2 

NO  * 2 H 20C0  2* 

N2 

2.00E-09  0.0 

50 

NO*2H20C02*  M20 

= 

NO* • 3H20 

4 

C02 

2.00E-09  0.0 

51 

N0*.N2 

♦ 

N2 

2 

NO* 

4 

N2 

4 

N 2 

1.00E-08  -4.  4-2. 10E*  03 

52 

H904* 

♦ 

H 

2 

H703* 

4 

H20 

4 

N 

1.20E-01  -4. 0-0.80E* 03 
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Eight  negative  ions  arc  included  in  this  program,  O ",  O ",  CO.  , CO„  , O , 
03",  NOg"  and  NOg".  This  is  also  the  order  in  which  the  respective  ion  concentra- 
tions are  solved.  Considerable  algebraic  manipulation  was  required  for  these  solu- 
tions which  will  not  be  described  in  the  text  but  can  be  followed  in  the  subroutine 
NEGJON,  lines  25-87.  The  ion  concentrations  of  several  ions  are  easily  solved,  but 
the  main  trick  was  to  eliminate  COg  from  the  equation  for  Og  which  was  accom- 
plished by  reducing  COg  to  a function  of  Og  . Thus,  the  solution  of  COg  (line  64) 
is  expressed  simply  as  S [Og  ] / R,  but  S,  lines  35-39  and  R,  line  34,  are  quite 
complex  terms  involving  the  lengthy  algebraic  factors  T,  W,  V,  U,  which  are 
described  on  lines  29-32.  The  negative  ion  processes  used  in  the  calculations  are 
given  in  Table  2.  Rate  coefficients  for  non-proton  processes  were  taken  mainly 
from  the  work  of  the  Boulder  group  (for  example,  Ferguson)10  including  a recent 
update11  for  kon  + k01,  8.  8X  10  10  cm3 /sec  with  product  branches  which  we  took 

oU  o 1 

as  equal.  Photodissociation/photodetachment  rates  are  generally  those  listed  by 

12  13- 

Peterson  et  al  or  estimated  by  Niles  (NO^  ). 

Profiles  of  neutrals  O,  Og,  Og,  Og(  1 A ' , Ng,  COg,  HgO,  NO  and  NOg  plus 
temperature  were  required  inputs  to  the  computations.  Total  concentrations,  [M], 
and  temperatures  T,  were  taken  from  CIRA  197  2 for  the  appropriate  month  and 
latitude  with  N0=  0.  7808  [M],  O = 0.  2095  [M]  and  CO„=  3X  10"’  [M].  Water  vapor 

concentrations,  [H0O]  were  set  at  10  [M].  The  choices  for  the  other  minor 

3 

neutrals  are  described  elsewhere.  Each  output  gives  a complete  account  of  all  the 
neutral  concentrations  utilized.  The  daytime  and  nighttime  minor  neutral  concen- 
trations used  in  the  computations  will  be  presented  in  typical  printouts  for  daytime 
and  nighttime  situations. 


10.  Ferguson,  E.  E.  (1975)  in  Atmospheres  of  Earth  and  the  Planets  (Ed.  B.  M. 

McCormac)  Reidel,  Dordrecht,  Holland,  pp  197-210. 

11.  Ferguson,  E.  E.  (1977)  Paper  presented  at  IAGA  meeting  in 

Seattle,  Washington, 

12.  Peterson,  J.  R.  , Cosby,  P.  C. , and  Moseley,  J.  T.  (1977)  Photo-destruction  of 

atmospheric  negative  and  positive  ions.  Space  Res.  F?:243-246. 

13.  Niles,  F.  E.  ( 1976)  private  communication. 

COSPAR  International  Reference  Atmosphere. 


Table  2.  Negative  Ion  Reactions.  Units  and  format  are  as  for  Table  1 
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3.  NUMERICAL  EXAMPLE:  2 November  1969  SPE 

Table  3 is  a typical  printout.  It  was  generated  specifically  for  the  ionization 
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production  rate  associated  with  the  daytime  electron  SPE  concentrations  measured 
on  2 November  1969,  Flight  B4  at  Ft.  Churchill,  Canada.  Photodetachment/ 
photodissociation  rates  kgl  - kg?  have  been  doubled  to  improve  the  agreement  be- 
tween the  data  and  the  computations.  The  agreement  is  quite  good,  as  seen  in 
Figure  1.  The  results  obtained  with  the  "normal"  Table  2 values  for  kgi  - kg?  are 
discussed  elsewhere3  in  conjunction  with  the  doubled  rates.  The  rate  coefficient 

k,_  has  been  set  to  0 in  this  example  because  NO+  otherwise  seemed  to  fall  off  faster 
43  15  + 

with  decreasing  altitude  than  the  data  for  this  event  admit.  This  NO  behavior 

is  shown  in  Figure  2 which  is  a plot  of  the  principal  positive  ions  listed  in  Table  3. 


14.  Dean,  W.A.  (1972)  in  Proc.  COSPAR  Symp.  on  SPE  of  Nov.  1969,  AFCRL-72- 

0474,  pp  291-306. 

15.  Narcisi,  R.  S. , Philbrick,  C.  R.  , Thomas,  D.  M. , Bailey,  A.  D.  , Wlodyka,  L.  E. , 

Baker,  D. , Federico,  G. , Wlodyka,  R. , and  Gardner,  M.  E.  (1972)  in 
Proc.  COSPAR  Symp.  on  SPE  of  Nov.  1969,  AFCRL-72-0474,  pp  421-432. 


Also  shown  in  Figure  2 is  an  NO+  profile  obtained  with  k^^  0,  specifically,  the 
the  value  given  in  Table  1.  The  output  associated  with  k^^  0 is  given  in  Table  4 
minus  the  first  part  of  the  output  which  is  the  same  as  that  for  Table  3.  The 
electron  concentrations  are  somewhat  higher  for  the  first  case,  the  greatest  differ- 
5 ence  being  near  7 8-80  km.  This  is  because  with  k ,„4  0,  O0  + /NO+  is  higher  at 

80  km  and  Og  has  a lower  dissociative  recombination  coefficient,  twice  that  of  NO 

(as  adopted  here.  Table  1).  In  other  words-  for  the  same  ionization  production 

. 1 / 2 + 

rate  the  electron  concentration  will  be  a factor  (2)  greater  if  only  09  ions  are 

present  as  compared  to  solely  NO  ions.  Slightly  lower  electron  concentrations 

are  encountered  near  75  km  for  the  k^^  0 case  because  conversion  of  NO+  to 

oxonium  ions  leads  to  slightly  higher  electron  loss  since  the  latter  ions  have  greater 

dissociative  recombination  coefficients  than  NO+  (Table  1).  Below  75  km,  there  is 

virtually  no  difference  in  [e]  for  the  two  cases,  because  NO  contributes  very  little 

to  the  ion  chemistry  in  either  case.  In  fact,  since  O^  is  the  precursor  ion  for  the 

SPE  situation,  the  choice  of  k^  is  rather  moot.  Furthermore,  inclusion  of  colli- 

sional  breakup  expressions  for  ions  like  NO+  (NgMCOg),  would  probably  increase 

the  NO+  concentrations  for  the  k^.^4  0 case  and  further  cloud  the  issue  as  to  the  role 

of  NO+  clustering  in  the  D- region  under  SPE  conditions.  Another  important  factor 

not  to  be  forgotten  in  regards  to  the  NO+  distribution  is  the  fact  that  it  depends 

strongly  on  the  NO  profile  chosen.  The  NO  profile  selected  is  based  in  part  upon 

2 

its  compatibility  with  the  data  and  is  probably  the  strongest  reason  why  our  results 
should  not  be  construed  as  giving  a major  clue  as  to  the  role  of  NO+  clusters  in 
the  D- region. 

7 

Vanderhoff  and  Niles  report  that  the  photodissociation  of  oxonium  is  negligible 
considering  the  small  cross-sections  for  the  process.  On  the  other  hand,  photo- 
dissociation cross-sections  for  O^4,  Og^HgO),  and  Og^COg)  are  large.  ^ We 
have  included  the  recommended  photodissociation  rates  for  the  first  two  species; 
processes  4 and  5 in  Table  1.  The  inclusion  of  these  rates  had  the  effect  of  lowering 
the  crossover  point  between  the  simple  ions  (NO+,  02+)  and  the  oxonium  ions  by 
about  2 km.  Thus,  the  effect  is  not  a major  one.  A negligible  further  change  in 
crossover  point  altitude  was  found  when  these  rates  were  doubled  to  account  roughly 
for  scattered  and  reflected  light  as  has  been  done  for  the  rates  kg^  - kg^.  Photo- 
dissociation of  Og^COg)  is  not  significant  in  the  D- region  because  this  species 
is  not  produced  effectively  at  these  altitudes. 


16.  Moseley,  J.  T.  , Cosby,  P.  C.  , and  Peterson,  J.  H.  ( 1977)  Laboratory  measure- 
ments of  photodissociation  and  photodetachment  cross-sections  of  atmospheric 
ions,  EOS  58:698. 


Although  the  daytime  electron  concentrations  for  3 and  4 November  1969,  be- 

3 

sides  the  2 November  1969  case  shown  here  agree  satisfactorily  with  the  data,  it 
should  not  be  assumed  that  the  negative  ion  problem  is  settled.  The  hydration  of 
negative  ions  has  been  ignored  (Figure  3),  and  this  is  undoubtedly  unrealistic.  The 
agreement  noted  therefore,  only  implies  that  the  negative  ion  chemistry  is  satis- 
factory in  regards  to  the  electron  distribution  in  daytime.  Thus,  negative  ions  like 
NO  ” and  NO.”  that  appear  to  have  negligible  photodetachment  rates,  may  be 

Ct  O 

hydrated  to  an  important  extent  but  have  no  impact  upon  the  electron  chemistry 

because  electrons  are  presumably  even  less  readilv  detached  from  NO  hydrates. 

3 x 

Formation  of  CO„  (H„0)  has  been  shown  to  be  unimportant  in  the  daytime  because 

<3  “ 2 

the  photodestruction  rate  is  large.  Even  sunset  does  not  appear  to  be  influenced 
by  the  omission  of  hydration  processes,  but  sunrise  is  another  matter.  The  major 
reason  for  a poor  match  of  the  November  1969  SPE  sunrise  data  with  time-dependent 
calculations  using  the  same  chemistry  as  here  is  quite  likely  due  to  the  considerable 
hydration  of  negative  ions  that  occurs  unhindered  by  photodissociation  processes 
at  night. 


Figure  1.  Measured  and  Computed  Electron 
Concentrations  for  Rocket  B-4  Launched  at 
2110  UT,  2 November  1969.  The  solar  zenith 
angle  was  83.  2°.  The  computed  electron 
concentrations  are  from  Table  4 
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Table  3.  Complete  Printout  of  Calculations  for  2110  UT.  2 November  1969  (Cont) 


Table  3.  Complete  Printout  of  Calculations  for  2110  UT,  2 November  1969  (Cont) 


T 
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Table  4.  Same  as  for  Table  3 Except  k,  „ as  in  Table  1 (Cont) 
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4.  NUMERICAL  EXAMPLE:  Nighttime 


Table  5 contains  charged  particle  concentrations  appropriate  to  SPE  conditions 

during  the  night  of  3 November  1969  at  Ft.  Churchill,  Canada.  The  ionization 

production  rates,  principal  constituents,  and  temperatures  are  identical  to  those 

used  in  our  report  dealing  with  a model  for  the  nighttime  steady-state  disturbed 

D-region.  Our  current,  more  advanced  model  results  are  not  extremely  different 

in  regards  to  [e]  as  compared  to  our  previous  effort.  The  major  difference,  a 

somewhat  lower  [e]  now  in  the  80-65  km  region  is  actually  in  better  accord  with  the 

data.  The  results  in  Table  5 were  obtained  with  k^:/  0.  Uncertainties  in  the  NO+ 

chemistry  and  the  chosen  NO  profile  notwithstanding,  it  still  is  worthwhile  to  point 

out  that  the  calculations  accomplished  with  the  Table  1 value  for  k,  „ resulted  in  an 
+ 2-3 

NO  concentration  at  7 0 km  of  2X  10  cm  , far  below  the  value  in  Table  5 which  is 

2 

in  good  agreement  with  the  data.  The  major  positive  ions  are  plotted  in  Figure  4. 
Below  about  60-65  km,  the  positive  ions  are  in  thermodynamic  equilibrium. 


Figure  4.  Principal  Positive  Ions  Computed  (Table  5) 
for  0605  UT,  3 November  1969,  the  Time  of  Launch 
for  Black  Brant  Rocket  17.758 
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The  negative  ion  populations  (Figure  5)  are,  undoubtedly,  as  noted  earlier, 
subject  to  considerable  hydration  at  night.  This  possibility  has  little  impact  on  [e] 
unless  the  hydrated  negative  ions  have  much  different  ion-ion  recombination  rates 
with  oxonium  ions  than  do  the  parent,  unhydrated  ions.  The  little  experimental 
evidence  in  this  regard  does  not  favor  this  possibility  which  is  attested  to  by  the 
rather  good  agreement  attained  between  model  and  data  for  [e]. 


I0Z  I03  I04  10* 

NEGATIVE  ION  CONCENTRATION  (cm"3) 


Figure  5.  Principal  Negative  Ions  Computed  (Table  5) 
for  0605  UT,  3 November  1969 


Table  5.  Complete  Printout  of  Calculations  for  0605  UT,  3 November  1969.  See  Table  3 caption  for  further 
description  of  the  table.  Rate  coefficient  k.7  was  set  to  zero 
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Table  5.  Complete  Printout  of  Calculations  for  0605  UT,  3 November  1969  (Cont) 
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5.  NUMERICAL  EXAMPLE:  August  1972  SPE 


The  final  example  presented  (Table  6)  is  for  one  instance  during  the  3 to  11 

August  1972  SPE.  The  temperature,  T,  and  total  molecular  concentration,  [M]  are 

for  1 August,  65°  N,  being  interpolated  from  the  C1RA  1972  models  for  60°  N and 

70°  N on  1 August.  The  ionization  production  rates  are  those  for  1508  UT  on 
17 

4 August.  Computed  [e]  agrees  very  well  with  data  taken  by  the  Chatanika  back- 

scatter  radar  at  62  km.  However,  above  this  altitude  the  model  predicts  electron 

17 

concentrations  2 to  3 times  greater  than  those  reported.  The  situation  is  more 
serious  below  62  km  where  model  electron  concentrations  fall  off  more  rapidly 

with  decreasing  altitude  than  do  the  data,  the  lowest  altitude  (50  km)  measurement 

4-3  3-3 

being  1 X 10  cm  whereas  we  compute  2.  5 X 10  cm  (Table  6). 

The  disagreement  is  disappointing  in  view  of  the  good  results  obtained  with  the 

model  and  the  2 to  5 November  1969  SPE.  However,  different  seasonal  factors 

may  influence  the  important  minor  neutral  species  in  such  a way  as  to  account  for 

the  differences.  Note  that  we  have  used  the  same  minor  neutral  concentrations  O, 

O.^,  NO,  NO2  and  02  (A)  as  previously,  but  N2>  02>  C02  and  H20  have  changed 

since  they  are  multiples  of  [M]  which  has  been  altered.  We  have  suggested^  that 

higher  H20  concentrations,  say  10  ppmv,  near  the  80  km  region  would  improve  the 

model/data  agreement  considerably.  The  rationale  for  higher  HgO  is  that  noctilucent- 

clouds  are  prevalent  in  the  high  latitude  summer  mesosphere.  The  situation  near 

50  km  is  much  more  serious  since  the  discrepancy  is  nearly  an  order  of  magnitude. 

A higher  0„,  and  hence  O concentration  would  help  alleviate  the  problem,  but  in 

4 19 

fact  there  is  evidence  that  may  have  declined  during  this  event  due  to  the  en- 
hanced formation  of  NO  resulting  from  proton  ionization  processes.  A decrease  in 
Og,  and  therefore  O,  not  to  mention  an  increase  of  NO,  and  hence  NOg  , NO^  , 
would  tend  to  make  matters  worse,  if  anything.  The  precision  of  the  electron  con- 
centrations measured  by  the  Chatanika  radar  below  60  km  should  be  considered  also. 

It  may  be  that  the  (e]  values  have  an  important  systematic  error  with  decreasing 
altitude  below  60  km.  We  note  that  no  data  is  given  below  50  km  indicating  that  this 
is  the  last  acceptable  data  point.  Deterioration  of  the  (e)-measurement  may  be 
setting  in  at  a somewhat  higher  altitude.  The  Chatanika  radar  has  been  able  to  make 
D-region  measurements  only  during  this  rather  intense  SPE. 

17.  Reagan,  J.  B.  , and  Watt,  T.  M.  (1976)  Simultaneous  satellite  and  radar  studies 

of  the  D-region  ionosphere  during  the  intense  solar  particle  events  of 
August  1972,  J.  Geophys.  Res.  8^1:4579-4596. 

18.  Swider,  W.  , and  Foley,  C.  I.  (1977)  The  influence  of  minor  atmospheric  con- 

stituents on  the  electron  loss  rates  of  the  August  4-11,  197  2 and  November  2-5 
1969,  solar  proton  events,  EOS,  58:465. 

wv» 

19.  Heath,  D.  F. , Krueger,  A.  J. , and  Crutzen,  P.  J.  ( 197  1 ) Solar  proton  event: 

Influence  on  stratospheric  ozone.  Science  197:886-889. 
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Reagan  and  Watt  claimed  that  xl/  (or  SI,  see  Eq.  (1))  for  this  event  agreed 
4 

with  the  results  for  the  2 to  5 November  1969  SPE.  However,  the  events  were 

not  normalized  in  terms  of  their  appropriate  T and  [M]  values.  The  l//  curves  for 

the  two  events  thus  agree  well,  save  for  the  80-km  region  problem,  only  if  T and 

[M]  for  the  August  197  2 SPE  was  similar  to  that  for  the  2 to  5 November  1969  event. 

This  possibility  appears  to  be  remote.  The  1 November  60°  N CIRA  1972  T and 

20 

[M]  values  used  for  the  November  SPE  are  corroborated  by  experimental  evidence. 

20.  Faire,  A.  C.  , and  Murphy,  E.A.  ( 197  2)  in  Proc.  COSPAR  Symp.  on  SPE  of 
Nov.  1969,  AFCRL-72-0474,  pp  445-455. 


Table  6.  Complete  Printout  of  Calculations  for  1508  UT,  4 August  19A2.  See  Table  3 caption  for  further 
description  of  the  table.  Rate  coefficient  k.,  was  set  to  zero.  The  solar  zenith  angle  was  7 9.  5° 
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Table  6.  Complete  Printout  of  Calculations  for  1508  UT,  4 August  1972  (Cont) 
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Appendix  A 


Steady-State  Multi-Ion  Disturbed  D- Region  Model  Program 
(Subroutine  Initial,  which  contains  input  profiles  of 
Q,  M,  T,  minor  neutrals  and  SI,  has  been  omitted.) 
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I 


40  - 90  KM  III  INCREMENTS  OF  2 


5 C 
C 

C 78  REACTIONS 

C 

PROGRAM  OREGION(INPUT=65,OIITPUT=65,TAPE15*65,TAPE2> 

10  COMMON  NEUT(9,26),TEMP(26I  ,M(26) ,Z(26) ,0(26)  ,SI (26) 

DIMENSION  NION(8 , 26)  ,NSPECI(8) ,PION(  30,26) ,<(100) ,ESAVE(26) 

1 , ISPECI  ( 30  > , INE'JT(9)  , A(100),  B <1 3 0 ) , 0(100), IREAC(IOO), 

2 JOUR ( 26) , I TER ( 26) , REAC (6) 

DIMENSION  ICOEFF (26)  ,EINT( 26) , SPI NT( 26) , ICO.  ( 2)  .EVENT (2) 

15  REAL  N, K,NEUT , NION 

INTEGER  HATER, TPF , PC  A 
REAL'  ICOEFF 

LOGICAL  IONCF, OAT, FIRST, OIAGNOS 
DATA'  ICOL/lOH  INT.  E ,10M  INT.  SP  / 

20  C 

C NALT*NUMBER  OF  ALTITUOES,  ISTART  IS  INDEX  1r  STARTING  ALTITUDE, 

C NELE=NUMBER  OF  ELEMENTS  IN  PION  ARRAY  TO  BE  COMPUTED 

NELE  = 3 0 

25  FIRST*. TRUE. 

C READ)  IN  COMPUTER  RUN  OPTIONS 

C 

READ  72,  ISTART, MALT 
30  72  FORNWM2I5) 

IF(EOF(5LINPtlT))  99,73 
73  READ  69, SET ,ISET 

69  FORMAT ( A 1 0 » 15) 

PRINT  69,  SET 

35  74  READ  70 ,OAV, MATER, TPF, PCA,  I ONCF, OIAG NOS.CRIT  , 00 ,E VENT ( 1) .EVENT (2) 

70  FORMAT ( 1L2, 3A2 ,2L2,2F5.0,2A10) 

IF(FOF (5LINPUT))  99,71 

71  PRINT  67, OAV.HATER, TPF, »CA, IONCF, 00,  GRIT, EVENTfll  ,EVENT<2) 

67  FORMAT ( * DAY=*,L2,»  H20**,A2,*  TPF=*,A2,*  PCA=*,A2,*  IONCF**, L2, 

40  1 * 00**»F5.0*  CRITERIA**, F5 .0'  PERCENT  * , 

2 7*  * , 2A10) 

C DETERMINE  SPECIE  PROFILES  AND  OTHER  PARAMETERS  WHICH  DECRI9E  THIS 

C COMPUTER  RUN  - FOUMO  IN  SU3R0UTINE  INITIAL 

45  C 

CALLI  INITIALIOAY, HATER, TPF, PCA) 

C 

PRINT  75  ,Z(ISTART),Z(NALT| 

75  FORMAT  < 5X ,* ALTITUDE  RANGE  *,  F5.0,*  -*,F5.),*  KM*) 

50  00*0. 

IF(DAY)  00*1. 

IVESM  S NOS  = MALT-ISTART*l 

* 

C 00*  OPTICAL  DEPTH  FACTOR 

55  C 

C COMPUTE  THE  IONIZATION  COEFFICIENTS  FOR  EACH  ALTITUOE 

P4 


SIGMA* 1 . 0E- 20 
IF(IONCF)  7,9 


60 


65 


70 


75 


SO 


85 


90 


95 


7 DO  6 J=ISTART,NALT 

H=(3032.*TEMP(J>>/(( 6370. '(6370.*  ZIJ)  )>**2) 

ICOEFF( J)=4.7E-07*E*P(-(NEUT(2, J> *H*SI GMA*03 1 ) 

PION(30, J)=ICOEFF( J)*NEUT(5, J) 

8 CONTINUE 
GO  TO  15 

(NO>*I  MILL  NOT  BE  3 ONUTEO  THIS  RUN 

9 NELF*29 

15  CONTINUE 
TOOAiY=10H 
RUNDY=DATE (TODAY) 

WRITE( ? , 150  I RUNDY.SET 

150  F0RHAT(1H1,15X,*RUN  OATE»  , » 10, 15X ,*RE ACTI ON  SET  0F*,A10) 
ALPHAO=4.0E-07 

ALPHAI=6.0E-08 

TIMF=10HOAY 

IF(  .NOT.  OAT)  TIME* 10HNTJHT 

MRITE(2, 151)  ALPHAO, ALPHAI, TINE, M ATER , TPF , P3 4 . CRIT, 

1 EVENT (1),EV'NT(2) 

151  F0RNAT(5X,*AD=*,1PES.1,»  A I =* , 1PE 8.1 , 2X, A10,  » , H20=*,A2, 

1 *,  TENP=*.A2,*,  ®CA=*,A2,*  CRITERIA=*,0PF5.0,*  PERCENT*, 

2 *,  EVENT  »,2A10) 

CRIT  =CPIT/ 1 03 • 

IF(  I'ONCF)  11,13 

11  MRI TE( 2 , 12)  30, SI  SNA 

12  FORMAT ( 10X,*3PTICAL  DEPTH  F ACTOR  = *,F5.2,*  SIGMA  = *,1PE9.1) 

13  NRFAC=78 

IF(  .NOT.  FIRST)  50  TO  19 

IF ( PC A .EQ.  2HL0  .OR.  »3»  .EQ.  2HHI  ) GO  T)  18 
IF( 0 CA  .EO.  2H72)  36,37 

36  RE  A IT  16,  (Q  ( I)  , 1=1 , 25  ) 

17  FORMAT ( 12F6.C ) 

16  FORN»T( 1P7E10.2) 

IF(EOF (5LINPUT) ) 99,18 

37  RE  A O'  17,  (Q(I)  , 1=1, 2>  ) 

IF(E  OF (5LINPUT) )99,18 

18  CONTINUE 


C REAT  IN  CHEMICAL  REACTIONS  FROM  T APE 1 5 


100 


105 


110 


00  20  1 1 = 1,  NR E AC 

REAO( 1 5,1)  I, (REAC(L),L=1,6),A(I) ,B(I) ,0(1) 

1 FORMAT (I4,A10,2A>,1X,A10,2A6, 1PE9  .2, 0PF5. l,l'>E9.2) 
IREAC(TI)=I 

IF(  = OF ( 15) ) 22,21 

21  IF( T SET  .EQ.  0)  50  TO  20 

MRI T E ( 2 , 2)  IREAC(t) , ( REAC ( L) ,L«1, 6) , A (I) , B(I ) , 5 (I) 

2 FORMAT (I4,A1Q,2A6,1H=,A10,2A6, 1PE9. 2 , 0PF5. 1,  1PE9. 2) 
20  CONTINUE 

GO  TO  19 

22  NREA3*I-1 
19  CONTINUE 

MRITE( 2, 805)  NREAC 

805  FOR’HT  ( • THE  NUMBER  OF  REACTIONS  THIS  RUN  *,110) 
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115  C IS®FCI  IS  A HOLLERITH  ARRAY  CONTAINING  CHEMICAL  SYM80L « OF  ®ION 

C ARRAlY  C POSITIVE  ION  ARRAT> 


C NSPrCI  IS  IOLLERITH  ARRAY  CONTAINING  CHEMICA-  SYMBOLS  OF 

120  C NIUN  ARRAY  (NEGATIVE  ION  ARRAY) 


IS®FCI  ( 1>=10HE 
ISPFCI < 2)  = 10HSP 
ISPECI t 3 ) = 1 0HNSU1 
125  ISPECI (4)=10H02» 

ISPECI (5)=10M04» 
ISPECI(6)=10HO2».H23 
ISPECI ( 7)=10HH3O» • ON 
ISPECI ( 8 )=10HH  30 Y 
130  ISPECI ( 9) * 10 HH30* .N 2 

ISPECI (10)  = 10HH  30* .C  02 
IS®ECI(11)= 10HH532* 
ISPFCK  12)  =10  MH50*.C02 
ISPECI ( 13) =10HNO» 

135  ISPECI  ( 14)  = 10HNO*  . N2 

ISPECI ( 15)  = IOHNOy . C3  2 
ISPECI  ( 1 6)  = 10  HNO* . H2  0 
IS®ECI(17)=10HNO+.H2O.N2 
ISPECI (1 8) =10NNO»H2O.CO2 
140  ISPECI ( 19) >10HNO».2N 20 

ISPECI  ( 20)  = 10HNO»  . 21 20N2 
ISPECI ( 21) = 1BHN012H20C02 
ISPECI (22)  = 10HNO* . 3M  20 
ISPECI (23)* 10HH7D3* 

145  IS®ECI ( 24)  = 10  HH934Y 

ISPECI ( 25) *10H  LA 
ISPECI ( 2 6)  = 10HL  AN30A 
ISPECI (27)=10HQ/E**2 
ISPECI  ( 28)  = 13MPStlM 
150  ISPECI!  29)  = 10HNSl)MlE 

ISPECI(30)=10H(NO) .1 


NS®ECI (1) =10HO2- 
NSPECt (2) =10H04- 
NS®ECI (3)=10HC04- 
NSPEd  ( 41  =10H03- 
NSPECt (5)=10HCO3- 
NSPECt (6) =10HO- 
NSPECt (7)*10HNO2- 
NSPECt ( 8) *10HNO3- 


C INEITT  IS  A HOLLERITH  ARRAY  CONTAINING  CHEMICAL  SYMBOLS  OF  NEUTRAL 

C SPECIES 

155  C 

INEITT  ( 1)  = 10H9 
INEITT (2)  *10H02 
INERT  ( 3)  « 10H3  3 
INEITT  (4)  *10 HI 2 
160  INERT(5)*10HNO 

INEITT  (6)  =10HC02 
INEITT  ( 7)  = 10HH20 
INEITT  (8)  =10  HS  2(13) 

INEITT(9)  = 10HN02 

165  C 

IF(.MOT.  OIAGNOS)  SO  TO  48 

WRITE (2, 271)  (ISPECI(J), J=25,29) ,IC3L(2),IS»ECI (2), 

1 ICOL(l) 1 1 SPEC  1(1) 

40  CONTINUE 

170  C 

C 

C ITER  ARRAY  OF  NUMBER  OF  ITERATIONS  TO  COMPUTE  E 

C 

C JO  II  N't  ARRAY  OF  SMICHES  INDICATING  CRITERIA  SATISFIED 
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1 


00  25  J=  1ST  ART  » MIL  T 
ITER (J) =0 
25  J0UN(J)=3H  NO 

00  LiOOP  100  CALCULATES  ALL  THE  POSITIVE  I0NJ 


50  00  100  J=ISTART,NALT 

IF( JOHN ( J)  .EQ.  3NVES)  SO  TO  100 
ITEP(J)  = ITER(J)U 

CALLI  R»TECN(<,A,9,D,  IREAC,  TEHP, J,  NR"  AC) 

IF( ITE R ( J)  . GT.  1)  SO  TO  90 

PI3NK25, J)=NEUT(2, J) • (K( 61 ) *NEUT ( 2, J ) ♦ K(62)» 
IF (WAY)  51,52 

51  PION(l,J)=S<JRT(Q(J)/SI(J)t 
GO  TO  53 

5 2 PI0N(1,  J)=SQRT(  t PI  OH  (25,  J)  7 (2.*»LPHAD>  >**20 
1 -PI0N(25,J)/(2.*RLPNA0) 

5 3 PION (2,  J)=SQRT( ( Q( J> * ( 5.*NEUT( 1, J)/A LPHAO^NE 
1 (5.*9EUT(1,J)*5EUTC3,  J) 

EINT(  J)  =PI0N(1,J( 

SPINT(J)=PI0N(2,J> 

90  AISP=ALPHAI*PI0H(2, J> 

CALL  NEGI0N(DAY,NI0N,K,PI0M(25,J) ,PI0N(1, J), 
AISN=ALPHAI*PI0N(3,J) 

C0N1  = K(8)*NE9T(1,  J)»K(9)  MEUTC7,  JI*K  < 4 2)  *HO 
C0N2=K(7)*NE9T(2»J)*M(J) 

C0N4*NEUT (7,J)*(K( 10) ♦Kllll ) ♦AISN+M 30)*PI01 
CON5*DO*K{4)MOO*K(5)*K(9>»NEUTC7,J)/CON4) 
PI0N(4,J)=,  9*0  tJ  ) 7(AISN*K(27)*PI0N(1,J)  *C0N? 
1 NEUT(5,J)-(C0N2'  l K<  8 ) *N  EUT  ( 1 , J)  ♦<( 

1 CCON1«X(28)»PION(1,  J>»)  ) 

PI0NI5, J)=PI0N(4, J)*  (CON27CCONl»K(28)*PION(t 
PI0N(6, J)=PI0N(5,  J)*K(9) *1EUT(7tJ)/; 0N4 
PION( 7 , J) =K ( 10) »NEUT(7,J)*PI0N(G, J) 7 
1 (K(t>2)*NEUT(7,J)«-AISN»K(30)» 

PION  (8, J) -PION (6, J)*K(11)*NEUT(7, J) 7 

1 (K 1 31 1 »PI0N(1,J) ♦ K(18)M 

2 K(13)*NEUT(4»J)**2*AIS'( 
C0N1.=K  (14)*NEUT(&,J> 

C0N?=K(32)*PI0N(1,J>  ♦ AISN 

PIONC  9,  J)=(K(13)»NEUT(4,J)*NEUT(4,J)  *PION(» 
1 OON1*K{35)*PION(1,  J)*AIS( 

C0N3*K (15) *NEUT (7,  J) 

PION (10,  J>  = (C0N1*PI0N(  9,J)  )/(C0N3+<  (36)*PI0 
H20NlxNEUT<7, JI*H  (J) 

PI0*(11, J)*(PI0N(7,J)*K( 12) *NEUT(7,J)*K(18>* 
1 P'ION(  10  » J)  *C0N3 ) 7 

1 (KC32)*PION(l, J»»K(19>» H20M1+AI 

1 K(16)*NEUT(6, J)  *N(J») 

PIONC 1 2 , J) s (5(16) *NEUT(&, J)*PION( 11,  J)*M( J)) 
1 (<(37)*»ION(l,J)»AISN*K(17 

C0N1  = K(21)*NEUT(5,  J)*NEUT(4,J> 
C0N2*KC22>*NEUT(7,  J) 

S0N2=NEUT (4,J)*NEUT(4»J) 


NE  l)T(  4,  J)  ) 


( J) /ALPHAD) 

JT(3,  J)/ALPHAI)(  7 
) ) 


J»AISP,PI0M3»J)  ) 
)*!JD*K{4>*AISN 
(1,  J)  ♦0D*K(5) 


H<C6>  * 

42)*N(J)»C0N5)/ 


PI0N(1, J)) 

EUT (7  , J)  *N  ( J)  ♦ 


, J>>7 


N(1,J)*AISM 
PI0N(8, J)*H20M1» 


I *NEUT(7, J)) 


230  PI3N(13, J)*(.1*Q( Jl  ♦<(&>» 1EUT< 5, J)*PI0N(4,J)  *PION(3Q,J))/ 

1 <KCZ6)*PrO*Cl, J)»C0N1»AISN 

2 ♦<  ( 43)  *S>N2-U<43)*SQN2*NEUT<4,J)  •K(51>' 

3 (K(44)*NE0T<6r J I »K < 38 I » ® I ON <1 , Jd AISNHU51)  »NEUT<4, Jill) 
PIONd4,  J)  = K(4  3)*PION(13,  J)  *S QN2)/ 

236  1 ( K(44)»NEUT<6,  J)  ♦K<38)  *PION  (1,  J)  ♦AISW  < (51)  *NEUT  (4,  JH 

PI0N<15,  J)  =<C0N1*PI0N(13,JI  ♦<  (44)  *NEUT  (6.  J)  ► » I ON( 14,  J>  )/ 

1 <CON2*K(JB)*PION<1,  J) +AISNI 

CONl*KC23)*NEUT(7,JI  *NEUT(4,J) 

PI0N(16,J)=C0N2*PI0N<15,  Jl / (CONl*K(3 9) *PI ON<  1 , J)» AISN 
240  1 ♦K(45I*SQN2) 

PIDN<17, J)=(PI0N<16,  J)*K(45)*SQN2)/ 

1 (K(45I*NEUT(6,JI*KC40)'  ®IONd,  JItAISN) 

PIONHie.JI*  CPI0N(17,  J)*K(45)*NEUT  (6,  Jll/ 

1 <<(47)  *NEUT<7, J)*<  (4C)»PI3Nd,  JI  + AISN) 

245  C0N2  = <(24)*NEUT<7,  J)  *NE’JTU,  J) 

PI0N<19, J)= (C0N1*PI0NC15, Jl ♦PIONl 18, J) *K(47) »NEUT(7,JI)/ 

1 (K<40>*PION(1»J>  ♦CON2*AISN»<  <48)*SQN2) 

PI3*(20,  J>  = (PI0N(19,  J>  *K(  48)  *SQN2)/ 

1 <K(49I  »NEUT(6,  J)  *K  (41)  *PI  ON  1 1,  J)  ♦AISH) 

25  0 PION(21,J)= (PION<20,  J)*K(49)*NEUT<6,  Jl)/ 

1 CKC50I  *NEUT(7,J) ♦ Kdl 1 * 3 ION<  1 , J) ♦ AISNI 

PI0N<22, J)=(CON2*PION(19, J) ♦ PIOIM 21, J) *K(50I  *NEUT(7,J))/ 

1 <K<25)*NEUT(7, J)  ♦ K<4 1 1 *PION(  1 , JltAISN) 

G=K(341*PIONd,J)«-AISN4*(5  2)*N(J) 

255  PI0NI23,  J)  = <PION(ll,  JI*KU9)*H20M1*PION(22,  J)  *K (25) *NEUT ( 7, Jl  ♦ 

1 Kd7)*PION(12,J)*NEUT<7,J))/ 

1 (K(3TI*PION(l,JI  ♦<<20>*H2OHl*Cl-K(52l*N<Jl  /GMAISN) 

PI0N(2  4,J) = <PION<2  3, J) *K (2C I*  H20M1) / G 
PI0N(26,  J)=PION<3,  Jl  /PION(l,  J) 

260  PI0N<27,  J)  = Q(  J) /PI  ON ( 1, Jl ** 2 

PI0NI29, J)*PlON(S,  Jl  ♦PIONC l . J) 

100  CONTINUE 

: IF( JUMP  .EQ.  2)  GO  TO  126 

265  C CON»UTE  THE  SUM  OF  THE  POSITIVE  IONS 

C 

OO  120  J*ISTART,NALT 
IFC  ITFR(J)  .GE.  201  J0JN(J)*3MN0 
IF( J DUN ( J)  .EQ.  3HVES)  SO  10  121 
270  ESA/'E  ( J)  =PI  ONd  • J) 

PIONf 28 , J) * 0. 

00  121  1*4, 24 

121  PION(26,JI*PION(28,J)4PION(I,J) 

IF  ( A8S(1.-PION(29,J)/PION(28,  J»)  E.  CRITI  JOUN <J)*3HYES 
275  IF(  JOUNIJ)  .EQ.  3HTES)  GO  TO  119 

PION <2,  J)=.5»  (ESAVE(  Jl  ♦PION  ( 28 , J)  ♦PTON(3,  J)» 

118  PION Cl, J)>.5»(ES1VE(  Jl  »PI0N(28«J)/(l. ♦PION I?  6, Jl  I I 
GO  TO  120 

119  ITESsITESd 

280  NOS=NOS-1 

120  CONTINUE 

IF  <WI  AG  NOS  140,49 
40  J*ISTART 

WRITE (2, 269 1 Z < Jl  , <»ION< I, Jl , 1*25 ,29) ,SPINT(  Jl  , PION (2, J) ,EINT(J) , 
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285 


1 


290 


295 


300 


305 


310 


315 


320 


325 


PION(l,J),ITER(J)  ..IDUNC  J» 

IF( ITER ( Jl  .SE.  20)  125,50 
99  CONTINUE 

IF (ITER (1ST ART)  .Ed.  1)  S3  TO  122 
00  125  J*ISTART,NALT 
IF(JOUN(J)  .EQ.  3H  90)  SO  TO  50 
125  CONTINUE 

PRINT*,-  THE  SCORE  IS  ",t»ES,"  YES  ANO  “,N3S,  “NO- 
122  WRITE (2,200)  (INEJT(I), 1*1,9) 

200  FORMAT ( /•  ALTITUOE* , 3X,9(A10) , 2X,*N*, 9X»  *TE  HP*) 

DO  235  J*ISTART, NALT 

WRITE( 2, 290 ) Z( J) , < NEUT( I, J) , 1*1, 9) , H ( J) , TE1  > ( J) 

290  FORNATIOPFIO.  0, 11 ( 1< , 1PE9. 2) ) 

235  CONTINUE 

WRITE (2,273)  ISPECK 1) ,IS*ECI ( 31 , (NSPECICJ), 3*1,8) 

273  FORHAT ( //*  ALTITUOE  *, 2X, 11 (IX , A10)l 
00  272  J*ISTART,NALT 

WRIT£( 2,279)  2(J) , ® ION(l, J) ,PION (3, J) , (NI01 (I , J) ,1 *1, 8) 

279  FOR*AT(F10.0,11(2X,1PE9.2)) 

272  CONTINUE 
IL*3 

00  270  1*1, 2 

19*11*1 

IL*IM*10 

1F(I1  .ST.  29)  IL* 29 

WRITE( 2,273)  (ISPECI ( J) , J*IM,IL) 

00  276  J*ISTART, MALT 

WRITE (2, 279)  Z(  J)  , (» I ONI N, J) , N*IH  ,IL) 

276  CONTINUE 

270  CONTINUE 

WRITEt  2 , 271)  (IS>ECI(J), J*25,29), ICOL ( 2) , IS* ECI (2), 

1 I30L(1),IS>ECI(1) 

271  FORMAT  ( //•  ALTITUDE  *, 2X , 9 ( IX , A10  ) , 2X , *Q* , 7< , * ITER*) 

00  275  J*ISTART,NALT 

WRITE (2, 269)  2 ( J)  , ( » I ON< I , J) , 1*25,29)  ,SPINT(  J)  , PI  ON (2, J) ,Et NT ( J) , 
1 »ION(l,J),Q(J> ,ITER(J> 

269  FORMAT ( OPF 10*  0,10(2X,  1PE9.  2)  ,2X,I5) 

275  CONTINUE 

IFdTER(ISTART)  .EQ.  1)  S3  TO  50 
FIRST*. FALSE. 

GO  TO  79 
99  CONTINUE 
CALL'  EXIT 
END 
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SUBROUTINE  REGION 


1 

SUBROUTINE 

C 

c 

SUBROUTINE 

5 

r * 

r* 

NE 

3ATI V E I 

c 

INDEX 

NT'ON 

•j 

1 

02- 

c 

2 

04- 

10 

3 

C04- 

c 

4 

0 3- 

«# 

5 

C03- 

c 

6 

0- 

r 

7 

N02- 

15 

8 

N03- 

20 


25 


S3 


35 


40 


45 


50 


55 


COMMON 
OINENST 
REAL  NT 
LOSI'CAL 
IFirMYJ 
0=1 

GO  TO  3 
0=0 

02SQ=NE 
X=X<  77) 
8=NEJT < 
1 <( 661 
T=tC(73» 
W=<  ( 71 ) 
V=K(68) 
U=<( 78) 

1 

R=0*'IU* 

s=o?sa* 

1 

2 

3 

4 


NEUT  IT  , 261 

ON  N10N(8,26I  , K <1001 
ON, K, LA, NEUT 
DAY 
1.2 


UT(2»J)*NEUT(2,J) 

•O2SQ«-K<80)  *MEUT<1,  J) 

1 , J) *( K (631  »< (64) ) HC  (65 )*NEUT ( 8 , J) ♦ 

*02SQ»0*K<9il ♦K<60>  * NEUT C3 , J) ♦ A ISP 

•NEUTC1,  J)»K<74)»NEUT<5, J) *K< 75 )*NEUT (9,  J) ♦0*K ( 54) *AI SR 
*NEUU1»J)*K(  72)*NEJT<6,J)  pAISP 

•NEUTI3,  J)»K<69>*NE'JT<5,  J)  »K<  70  )*NEUT  < 1, J) *0*K <93)  *AISP 

•NEUTIl, JI ♦ X*D*K( 9?  I *<(79) *NEU T <8, J). AISP 

♦K(8 1) *NEJT [3 1 Ji 

W*T-D»K<94I  *M72I • NEUT< 6, J )*X ) 

<W*U*X<66I*K<70>»NEJT<1, J) ♦ U’K  <66)*K< 58)*X<72I* 

NEUT <6, J)*NEJTC3, J)»0*K< 64) »K< 72 I •* (77)  *NEUT < 1, J) • 
NEUT (6, J)  ) 

♦V*K<72>  MEUTCd,  J)*NEUT<3,  J)  *C  J»Kl  50)  *K  <64)  *K<  80)  * 

N-  JT ( 1,  J) ) 


NIONd  , J)  =02- 

NIONC1, J)=< V*U*N*R*E*LA) /<R*9*V*U*N 

1 -J•R*K<65>*K<&8^ *<(71) *NEU T <1,  J) » MEUTC 3,  J)*02SQ- 

2 0*R*K<64)*K(71)*X*NEUTC1,  J)*»2- 

3 S*V*NEUT <l,J)*<K<71)*X*a*KC94)d*M*KC73)>- 

4 U*R*N*K<50)TK<71) • NEUT  (3 , J ) *NEUT  ( 1 , J) 

1 -H*V»K<64I  •K<81)  »NEUTC1,J>  *NEUT<3,  J) 

6 -<S*M*0*K<  81)*NEUT<3, J)*0*K<94) )/R> 

NI3N<2, J) =04- 

NION<2,  J)  = (K(66)  •02S<J*NI0N(1,J))7 
1 <K<67I*NE JT <6, JI  I 

NI0N<3,J)*C04- 

NION (3, J )* (Kt  66) *02SQ*  NI  ON  ( 1 , J ) ) / 

1 (K(68) *NEUT ( 3, J) ♦K<69) *NEUT<5, JI  *K< 70)  »NEUT<1 , J) 
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60 


2 


► D*K<931  ^ISPl 


C 

66  C 


7 0 


C 


nio»c;,j»*coj- 

NIO»<5, J)=S*NI0NI1,J)/R 
NI0N<6,J)=0- 

NI0«<6.  J>*  <K<64IMEJT<  1.JI  •NION(l.JMD<M(( 96)  » NT  ON  (5  , J>> /[) 
NI0<mfJI*03- 

NI0NI4, J)  = <K<661  »NEJT<3, Jl »NION<3,J) *K  <77 ) *41  OH (6 . J) *02SQ» 

1 NiUTIJ, Jl* <K<60IMION<1( Jl  * K < 60 1 » NION <6t Jl 1 1/N 


75  C 
C 


60 


C 


85 


90 


NI0N17, J1*N02- 

C=<0>*<  196)  )/(0*(K(96)  ♦ !(<  971  1 ♦ A ISP) 

NI0N<7,  J>*(KI7<*)*NEJT(5,  J)  »NI0N<5,  Jt  *C* « <75 1 *NE JT<9, Jl* NION < 5, Jl ♦ 

1 KI691* NEUT C5, Jl*NION<3, J1117 

2 <K<7S1MEUT<  3 » Jl  *0*K<951 ♦ AISP-C*K<75 1 *NEUT<3,  JH 

NIOMU, J1*N03- 

NIONJ6, Jl*<K<751*NEJT<9t J»  *NI0N<5,JI *K <761  *1EUT <3 ( J 1 *NION<7 , J 1 ♦ 

1 K< 69) ME9TI5.J1 *NION<3, Jll  • 

2 (0*<K<961*K<971 1.  *6ISPI 
SUN=.0. 

00  10  I*l«6 
10  SUN  “SUN*  NI  ONI  1 1 J) 

return 

END 


SUBROUTINE  RATECN 


1 SUBROUTINE  RA TECNI K, A , 8, 0, IREAC. TEN® , IALT , NR *AC 1 

DIMENSION  Kill fA<ll (Bill ,0<1I , IRE ACC  II , TEMPI  11 
REAL!  K 

00  10  J* 1(NRE AC 
5 I0*I«EAC<J1 

K<I  TO  ■ A < 101* ITEMP  (I  ALT! 73  00.1  **B  <101  *EXP  < 01  101  /TENPlIALTll 
10  CONTINUE 
RETURN 
ENO 
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DNA  DISTRIBUTION  LIST 


Department  of  Defense 
Director 

Defense  Advanced  Research  Projects  Agency 
1400  Wilson  Boulevard 
Arlington,  Virginia  22209 

1 cy  ATTN:  TIO 
1 cy  ATTN:  STO 
1 cy  ATTN:  NRMO 

Director 

Defense  Communications  Agency 
8th  Street  and  Courthouse  Road 
Arlington,  Virginia  22204 

3 cys  ATTN:  MEECN  Office 

Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virignia  22314 

12  cys  ATTN:  TC 

Director 

Defense  Nuclear  Agency 


Washington,  D.  C. 

20305 

1 cy  ATTN: 

STTL 

1 cy  ATTN: 

DOST 

3 cys  ATTN: 

RAAE 

1 cy  ATTN: 

RAEV 

Joint  Chiefs  of  Staff 
Department  of  Defense 

Washington,  D.  C. 

20301 

1 cy  ATTN: 

J-6 

Director 

National  Security  Agency 

Fort  George  G.  Meade,  Maryland  20755 

2 cys  ATTN:  Technical  Library 

Under  Secretary  of  Defense  (Research  and  Engineering) 
Department  of  Defense 
Washington,  D.  C.  20301 


2 cys  ATTN : 


ODS&SS 


Deoartment’of  Commerce 


U.  S.  Department  of  Commerce 

Office  of  Telecommunications 

Institute  for  Telecommunication  Sciences 

National  Telecommunications  and  Information  Administration 

Boulder,  Colorado  80303 

2 cys  ATTN:  W.  F.  Utlaut 


Department  of  the  Amy 

Command er/0 i rec tor 

Atmospheric  Sciences  Laboratory 

U.  S.  Army  Electronics  Conmand 

White  Sands  Missile  Range,  New  Mexico  88002 

1 cy  ATTN:  DRSEL-BL-SY-S 

F.  E.  Niles 

Director 

U.  S.  Army  Ballistic  Research  Laboratories 
Aberdeen  Proving  Grounds,  Maryland  21005 

1 cy  ATTN:  George  E.  Keller 

Commander 

U.  S.  Army  Foreign  Sciences  and  Technology  Center 
220  7th  Street,  N.E. 

Charlottesville,  Virginia  22901 

1 cy  ATTN:  Robert  Jones 


Department  of  the  Navy 

Chief  of  Naval  Operations 
Department  of  the  Navy 
Washington,  D.  C.  20350 

1 cy  ATTN:  NOP  985 

1 cy  ATTN:  MOP  094H 

Chief  of  Naval  Research 
Department  cf  the  Navy 
800  North  Quincy  Street 
Arlington,  Virginia  22217 

1 cy  ATTN:  Code  465,  R.  G.  Joiner 

1 cy  ATTN:  Code  427,  H.  Mu  Haney 
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> Comnander 

Naval  Electronic  Systems  Coumand 
Department  of  the  Navy 
Washington,  D.  C.  20360 


1 cv 

4TTN 

PME-117 

1 cy 

ATTN 

PME-117T 

1 cy 

ATTN 

PME-117-21 

1 cy 

ATTN 

PME-U7-21A 

1 cy 

ATTN 

PME-117-22 

Director 

0 

Naval  Ocean  Systems  Center 

Electromagnetic  Propagation  Division 

271  Catalina  Boulevard 

Sa n D.i ego , Cali f or n i a 92152 

1 cy 

ATTN 

Code  2200,  W.  F.  Moler 

1 cy 

ATTN 

Code  2200,  11  an  Rothmuller 

1 cy 

ATTN 

Code  2200,  John  Bickel 

Director 

Naval  Research  Laboratory 
4555  Overlook  Avenue,  S . W . 

Washington,  D.  C.  20375 

Code  7700,  Timothy  P.  Coffey 
Code  7709,  Wahab  Ali 
Code  7750,  John  Davis 
Code  2627 

Commander 

Naval  Surface  Weapons  Center  (White  Oak) 

Silver  Spring,  Maryland  20910 

1 cy  ATTN:  Technical  Library 

Office  of  Naval  Research  Branch  Office 
1030  East  Green  Street 
Pasadena,  California  91106 

* 1 cy 


cy 

1 cy  Ann 

2 cys  ATTN 
1 cv  ATTN 


Department  of  the  Air  Force 
Commander 

Air  Force  Geophysical  Laboratory,  AFSC 

L.  G.  Hanscom  Air  Force  Base,  Massachusetts  01731 


1 

cy 

ATTN: 

OPR, 

James  Ulwick 

1 

cy 

ATTN: 

LKB, 

W.  Swider 

1 

cy 

ATTN: 

L&B, 

K.  Champion 

r 


Director 

Air  Force  Technical  Applications  Center 
Patrick  Air  Force  Base,  Florida  32920 

1 cy  ATTN;  TD 

1 cy  ATTN:  HO  1035th  TCHOG/TFS 


Department  of  Defense  Contractors 

General  Electric  Company 

TEMPO  - Center  for  Advanced  Studies 

816  State  Street 

Santa  Barbara,  California  93102 

1 cy  ATTN:  Warren  S.'  Knapo 
1 cy  ATTN:  DAS  I AC 

Lockheed  Missiles  and  Space  Company 

3251  Hanover  Street 

Palo  Alto,  California  94304 


1 

cy 

ATTN: 

J.  B.  Reagan 

1 

cy 

ATTN: 

W.  Imhof 

1 

cy 

ATTN: 

Martin  Walt 

Mission  Research  Corporation 
735  State  Street 

Santa  Barbara,  California  93101 

1 cy  ATTN:  M.  Scheibe 

1 cy  ATT?!:  D.  Sowle 

Pacif ic-Sierra  Research  Corporation 
1456  Cloverfield  Boulevard 
Santa  Monica,  California  90404 

1 cy  ATTN:  E.  C.  Field 

Pennsylvania  State  University 
Ionospheric  Research  Laboratory 
College  of  Engineering 
318  Electrical  Engineering  - East  Wing 
University  Park,  Pennsylvania  16802 

1 cy  ATTN:  John  S.  Nisbet 

1 cy  ATTN:  Les  Hale 

1 cy  ATTN:  A.  J.  Ferraro 

1 cy  ATTN:  H.  S.  Lee 


44 


R&D  Associates 

4640  Admiralty  Way 

Marina  Del  Rey,  California  90291 

1 cy  ATTN:  R.  Lelevier 
1 cy  ATTN:  F.  Gilmore 
1 cy  ATTN:  R.  Turco 

The  Rand  Corporation 

1700  Main  Street 

Santa  Monica,  California  90406 

1 cy  ATTN:  Cullen  Crain 

Professor  Chalmers  F.  Sechrlst 
155  Electrical  Engineering  Building 
University  of  Illinois 
Urbana,  Illinois  61801 

1 cy  ATTN:  C.  Sechrist 

Stanford  Research  Institute 
333  Ravenswood  Avenue 
Menlo  Park,  California  94025 

1 cy  ATTN:  Allen  M.  Peterson 
1 cy  ATTN:  Ray  L.  Leadabrand 


